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We present a novel extrapolation scheme for high order series expan- 
sions. The idea is to express the series, obtained in orders of an external 
variable, in terms of an internal parameter of the system. Here we ap- 
ply this method to the 1-triplet dispersion in an antiferromagnetic S = 
Heisenberg ladder. By the use of the internal parameter the accuracy of 
the truncated series is enhanced tremendously. 

PACS numbers: 75.40.Gb, 74.25.Ha, 75.10.Jm, 02.30.Mv 

1. Introduction 

High order series expansions have become a powerful tool in the field 
of strongly correlated electron systems [1]. Especially in low dimensional 
quantum spin systems significant progress has been achieved, see e.g. [2-6]. 
Usually the obtained truncated series must be extrapolated using various 
extrapolation schemes in order to obtain results for the physically interest- 
ing regions [7]. Albeit these extrapolations are a standard technique they 
always introduce some uncertainty about the results. Therefore, a general 
transformation scheme which allows to read off the information content of 
the high order series more directly is highly desirable. Here we will pro- 
pose such a general scheme. Its usefulness will be demonstrated for the 
1-triplet dispersion in antiferromagnetic 2-leg S = ^ Heisenberg ladders. 
The series was obtained previously by various techniques [2, 4, 5]; we used 
a perturbative continuous unitary transformation (CUT) [3]. 

The Hamiltonian for the antiferromagnetic 2-leg Heisenberg ladder reads 



H 



E [-^11 (Sl,iSl,i+l + S2,iS2,i+l) + J±Si,iS2,i] (1) 



where Jy and J± are the exchange couplings on the legs and on the rungs, 
respectively. The subscript i denotes the rungs and the subscript 1,2 the 
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legs of the ladder. 

We use a CUT to map the Hamiltonian H to an effective Hamiltonian 
i?eff which conserves the number of rung-triplets, i.e. \Hq,Hcs\ = where 
Hq := i?|j||=o [3]. The ground state of i^eff is the rung-triplct vacuum. 
The effective Hamiltonian H^^ is calculated in the 1-triplet subspacc of the 
Hilbert space to order 14 in x := J\\/J±- The ground state energy Eq = 
(0|i7cff|0) and the 1-triplet dispersion u>{k) = {k\Hcf[\k) — Eq is obtained 
[2, 4, 5]. The 1-triplet dispersion is expressed in terms of the 1-triplet 
hopping amplitudes tnix) = (i|i?effK + where \i) denotes the state with 
one triplet on rung i 

oj(k)/ J± = ^'^tn{x) cos{nk) . (2) 

n 

The 1-triplet dispersion has a minimum for k = tt, i.e. the 1-triplet gap 
A(x) given by co{Tr). In Fig. la, the results for A{x) are shown. The 
truncated series for A{x) yields a satisfactory agreement up to a; « 0.6. 
Beyond this value it cannot be used as estimate for the value of the gap. 
More sophisticated extrapolation schemes, however, still work fine. 




Fig. 1. (a) Gap A{x)/J±. The solid and the dotted curves curves result from a 
[7,7], [6,8] and a [8,6] dlogPade approximant, rasp., assuming constant behavior on 
X ^ oo. The long-dashed curve depicts the truncated series. 

(b) Dispersion relation u{k) in units of J± + Jy for the couplings x = J\\/J± = 
0.2, 0.4, . . . , 1.6. The highest plot at fc = tt is a; = 0.2 and the lowest plot is 
X = 1.6. The long-dashed curve is the limit of isolated chains (a; = oo), i.e. the des 
Cloizeaux-Pearson result (7r/2) sin(A;) [10]. 
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2. Extrapolation Scheme and Results 

Generically, the various physical quantities in a given system depend 
in a complicated way on the external control parameters. Expanding the 
quantities under study in terms of one of the external control parameters, 
let us say yields the bare, truncated series which can only rarely be di- 
rectly used to compute the quantities. This is so since singularities induced 
by phase transitions easily spoil the convergence of the series. For instance, 
a correlation length diverges and the corresponding energy gap closes ren- 
dering an expansion about the gapped phase difficult. 

If the convergence of the series is deteriorated due to an incipient phase 
transition it is reasonable to assume that all quantities in the particular sys- 
tem behave in a similar fashion. If this is so we may proceed in two separate 
steps: (i) We extrapolate an internal parameter which may serve as a mea- 
sure of the distance to the phase transition as reliably as we can. Thereby, 
we can attribute reliably to a given value of x the corresponding value of the 
internal parameter, (ii) We express all other quantities as functions of the 
internal parameter. According to our argument the latter functional depen- 
dencies are expected to be much simpler, i.e. they are much less singular. 
The canonical candidate for the internal parameter measuring the distance 
to a phase transition or, more generally, to some singular situation is the 
energy gap A. It is inversely proportional to the correlation length ^ which 
plays the role of the internal control parameter in standard renormalization 
group treatments. 

In the following, we illustrate for the 1-triplet dispersion of the model 
in (1) that this scheme works indeed stunningly well. The gap A(a;) can 
be extrapolated reliably up to a; « 2 using dlogPade-approximants, Fig. la. 
In this extrapolation we can exploit additional properties of the gap such 
as its positivity and its asymptotic behavior for a; — oo. In this way, very 
reliable extrapolations are possible (Fig. la) so that step (i) is successfully 
solved. Note that A(x) in units of J_l does not vanish on x — > oo but rests 
finite [8, 9]. 

For step (ii) we define 

p{x) = 1 - A(x)/((1 + x)J^) = 1 - A(x)/(J|| + J^) . (3) 

In units of Jy -|- J_\_ the gap is unity at a; = and it goes to zero on a; — > oo. 
So p{x) varies monotonically between and 1 when x is increased from 
to oo. Note that p = 1, i.e. x = oc, constitutes the limit where the spin 
ladder becomes a system of two isolated spin chains. Since one has p oc x 
for small x any expansion in x can be rewritten as expansion in p of the 
same order as the series in x, yet with other coefficients! This is done by 
inverting Eq. 3, thus completing the second step. 
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We applied the extrapolation scheme proposed above to the expansions 
of the 1-triplet hopping amplitudes t„ in Eq. 2 which yields the 1-triplet 
dispersion Lj{k,p). In Fig. lb the results are depicted for the truncated 
series in p without any further extrapolation! For x < 0.6 the dispersion 
is a monotonic decreasing function in k whereas for larger values of x a 
characteristic dip at k = appears. For comparison, the limiting case of 
isolated chains is also included [10]. The main point is that the truncated 
series in p gives a quantitatively correct 1-triplet dispersion up to a; « 1.2 
(cf . Ref . [2] ) and qualitatively good results up to x 2 while the truncated 
series in x can be trusted only up to x 0.6. 

3. Conclusion 

We presented a novel generally applicable extrapolation scheme for high 
order series expansions. The basic idea is to use an internal control pa- 
rameter instead of an external control parameter as expansion variable. All 
difficulties stemming from singularities arc dealt with in the determination 
of the dependence of the internal parameter on the external one. 

The power of the novel scheme was demonstrated for the 1-triplet dis- 
persion of the antifcrromagnetic S = ^ 2-leg Hcisenberg ladder. Further 
investigations on other quantities and systems are in progress. In addition, 
the use of standard extrapolation techniques on the series in the internal 
parameter constitutes a very promising route to extend the applicability of 
strong coupling expansions. 

We are indebted to H. Monien and E. Miiller-Hartmann for helpful dis- 
cussions and to the DFG for financial support in SP 1073 and in SFB 608. 
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